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Abstract Surface parameterization is widely used in computer graphics and geometry processing. It 
simplifies challenging tasks such as surface registrations, morphing, remeshing and texture mapping. 
In this paper, we present an efficient algorithm for computing the disk conformal parameterization 
of simply-connected open surfaces. A double covering technique is used to turn a simply-connected 
open surface into a genus-0 closed surface, and then a fast algorithm for parameterization of genus-0 
closed surfaces can be applied. The symmetry of the double covered surface preserves the efficiency 
of the computation. A planar parameterization can then be obtained with the aid of a Mobius trans¬ 
formation and the stereographic projection. After that, a normalization step is applied to guarantee 
the circular boundary. Finally, we achieve a bijective disk conformal parameterization by a compo¬ 
sition of quasi-conformal mappings. Experimental results demonstrate a significant improvement in 
the computational time by over 60%. At the same time, our proposed method retains comparable ac¬ 
curacy, bijectivity and robustness when compared with the state-of-the-art approaches. Applications 
to texture mapping are presented for illustrating the effectiveness of our proposed algorithm. 

Keywords Disk conformal parameterization • Simply-connected open surface • Quasi-conformal 
theory 


1 Introduction 

With the advancement of computer technologies and 3D acquisition techniques, 3D objects nowa¬ 
days are usually captured and modeled by triangular meshes for further usages. A large number of 
applications of triangular meshes can be found in computer graphics and computer-aided design. 
However, working on general meshes is a difficult task because of their complicated geometry. The 
complicated geometry hinders applications such as surface registration, morphing and texture map¬ 
ping. To overcome this problem, one common approach is to parameterize the surfaces onto a simple 
parameter domain so as to simplify the computations. For instance, textures can be designed on the 
simple domain and then mapped back onto the original surfaces [ 1511231 157 ] . Another example that 
usually makes use of parameterization is surface registration hm. Instead of directly computing 
the registration between two convoluted surfaces, one can perform the registration on the simple 
parameter domain, which is much easier. It is also common to perform surface remeshing with the 
aid of parameterizations [ 7711251150 ] . With the development of the computer industry, the problem of 
finding a good parameterization method is becoming increasingly important. 
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To make a parameterization useful and applicable, one should seek for a method that minimizes 
certain types of distortions. In particular, it is desirable to minimize the angular distortions of the 3D 
meshes. Angle preserving parameterizations , also known as conformal parameterizations , effectively 
preserve the local geometry of the surfaces. Hence, in this paper, our goal is to develop an efficient 
conformal parameterization algorithm. 

The choice of the parameter domain is also a key factor in deciding the parameterization scheme. 
For simply-connected open surfaces, one popular choice of the parameter domain is the unit disk. 
Using the unit disk as a parameter domain is advantageous in the following two aspects. Firstly, the 
existence of the conformal parameterization is theoretically guaranteed. By the uniformization theo¬ 
rem, every simply-connected open surface is conformally equivalent to the open unit disk. Secondly, 
unlike free and irregular shapes on the plane, a consistent circular boundary of the parameter domain 
facilitates the comparisons and mappings between different surfaces. 

In real applications, besides the quality of the parameterization result, it is also important to 
consider the computational efficiency of the parameterization algorithm. In particular, a fast algo¬ 
rithm is desired so that the computation can be completed within a short time. In this paper, we 
develop an efficient algorithm for the disk conformal parameterization of simply-connected open sur¬ 
faces. To achieve the efficiency, we first transform a topologically disk-like surface to a genus-0 closed 
surface by a double covering technique. Then we can apply a fast parameterization algorithm for 
genus-0 closed surfaces to compute a spherical parameterization of the double covered surface. Note 
that although the size of the problem is doubled by double covering, the computational efficiency is 
preserved because of the symmetry of the double covered surface. The spherical parameterization, 
together with a suitable Mobius transformation and the stereographic projection, provides us with an 
almost circular planar parameterization for the original surface. A normalization technique followed 
by a composition of quasi-conformal maps are then used for obtaining a bijective disk conformal pa¬ 
rameterization. The bijectivity of the parameterization is supported by quasi-conformal theories. The 
entire algorithm only involves solving sparse linear systems and hence the computation of the disk 
conformal parameterization is greatly accelerated. Specifically, our proposed method speeds up the 
computation of disk conformal parameterizations by over 60% while attaining accuracy comparable 
to the state-of-the-art approaches mm- In addition, our proposed method demonstrates robustness 
to highly irregular triangulations. 

The rest of the paper is organized as follows. In Section[2j we review the previous works on surface 
parameterizations. In Section[3j we highlight the contribution of our work. Our proposed algorithm is 
then explained in details in Section [4j The numerical implementation of the algorithm is introduced 
in Section [5] In Section [6j we present numerical experiments to demonstrate the effectiveness of our 
proposed method. The paper is concluded in Section [7] 


2 Previous works 

With a large variety of real applications, surface parameterization has been extensively studied by 
different research groups. Readers are referred to [MUMS] for surveys of mesh parameterization 
methods. In this section, we give an overview of the works on conformal parameterization. 

A practical parameterization scheme should retain the original geometric information of a surface 
as far as possible. Ideally, the isometric parameterization, which preserves geometric distances, is the 
best parameterization in the sense of geometry preserving. However, isometric planar parameteri¬ 
zations only exist for surfaces with zero Gaussian curvature [5]. Hence, it is impossible to achieve 
isometric parameterizations for general surfaces. A similar yet far more practical substitute is the 
conformal parameterization. Conformal parameterizations are angle preserving, and hence the in¬ 
finitesimal shape is well retained. For this reason, numerous studies have been devoted to surface 
conformal parameterizations. 

The existing algorithms for the conformal parameterizations of disk-type surfaces can be divided 
into two groups, namely, the free-boundary methods and the fixed-boundary methods. For the free¬ 
boundary methods, the planar conformal parameterization results are with irregular shapes. In P’6], 
Hormann and Greiner proposed the MIPS algorithm for conformal parameterizations of topologi- 
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Methods 

Boundary 

Bijective? 

Iterative? 

Shape-preserving [6] 

Fixed 

Yes 

No 

mips [E! 

Free 

Yes 

Yes 

ABF/ABF++ [331133] 

Free 

Local (no flips) 

Yes 

LSCM/DNCP 31(23] 

Free 

No 

No 

Holomorphic 1-form [12] 

Fixed 

No 

No 

Mean-value 0 

Fixed 

Yes 

No 

Yamabe Riemann map [26j 

Fixed 

Yes 

Yes 

Circle patterns [22] 

Free 

Local (no flips) 

Yes 

Genus-0 surface conformal 

Free 

No 

Yes 

map [19] 

Discrete Ricci flow [20 

Fixed 

Yes 

Yes 

Spectral conformal [2? 

Free 

No 

No 

Generalized Ricci flow [36] 

Fixed 

Yes 

Yes 

Two-step iteration [3] 

Fixed 

Yes 

Yes 


Table 1 Several previous works on conformal parameterization of simply-connected open surfaces. 


cally disk-like surfaces. The boundary develops naturally with the algorithm. In [33], Sheffer and de 
Sturler proposed the Angle Based Flattening (ABF) method to compute conformal maps, based on 
minimizing the angular distortion in each face to the necessary and sufficient condition of a valid 
2D mesh. Sheffer et al. [33] extended the ABF method to ABF+-1-, a more efficient and robust 
algorithm for planar conformal parameterizations. A new numerical solution technique, a new re¬ 
construction scheme and a hierarchical technique are used to improve the performance. Levy et al. 
[23] proposed the Least-Square Conformal Maps (LSCM) to compute a conformal parameterization 
by approximating the Cauchy-Riemann equation using the least square method. In [3], Desbrun et 
al. proposed the Discrete, Natural Conformal Parameterization (DNCP) by computing the discrete 
Dirichlet energy. In [22], Kharevych et al. constructed a conformal parameterization based on circle 
patterns, which are arrangements of circles on every face with prescribed intersection angles. Jin et 
al. [19] applied a double covering technique m and an iterative scheme for genus-0 surface confor¬ 
mal mapping in [12] to obtain a planar conformal parameterization. In E2, Mullen et al. reported a 
spectral approach to discrete conformal parameterizations, which involves solving a sparse symmetric 
generalized eigenvalue problem. 

When compared with the free-boundary approaches, the fixed-boundary approaches are advanta¬ 
geous in guaranteeing a more regular and visually appealing silhouette. In particular, it is common 
to use the unit circle as the boundary for the conformal parameterizations of disk-type surfaces. Nu¬ 
merous researchers have proposed brilliant algorithms for disk conformal parameterizations. Floater 
[6] proposed the shape-preserving parameterization method for surface triangulations by solving lin¬ 
ear systems based on convex combinations. In [7], Floater improved the parameterization method 
using a generalization of barycentric coordinates. In na, Gu and Yau constructed a basis of holo- 
morphic 1-forms to compute conformal parameterizations. By integrating the holomorphic 1-forms 
on a mesh, a globally conformal parameterization can be obtained. In [26], Luo proposed the combi¬ 
natorial Yamabe flow on the space of all piecewise flat metrics associated with a triangulated surface 
for the parameterization. In [20], Jin et al. suggested the discrete Ricci flow method for conformal 
parameterizations, based on a variational framework and circle packing. Yang et al. [36] generalized 
the discrete Ricci flow and improved the computation by allowing two circles to intersect or separate 
from each other, unlike the conventional circle packing-based method [20]. In [3], Choi and Lui pre¬ 
sented a two-step iterative scheme to correct the conformality distortions at different regions of the 
unit disk. Table jT] compares several previous works on the conformal parameterizations of disk-type 
surfaces. 

Our proposed algorithm involves a step of spherical parameterization. Various spherical param¬ 
eterization algorithms have been developed in the recent few decades, such as 
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Among the existing algorithms, we apply the fast spherical conformal parameterization algorithm 
proposed in [2]. More details will be explained in Section [4] 


3 Contributions 

In this paper, we introduce a linear formulation for the disk conformal parameterization of simply- 
connected open surfaces. Unlike the conventional approaches, we first find an initial map via a pa¬ 
rameterization algorithm for genus-0 closed surfaces, with the aid of a double covering technique. The 
symmetry of the double covered surface helps retaining the low computational cost of the problem. 
After that, we normalize the boundary and apply quasi-conformal theories to ensure the bijectivity 
and conformality. Our proposed algorithm is advantageous in the following aspects: 

1. Our proposed method only involves solving a few sparse symmetric positive definite linear systems 
of equations. It further accelerates the computation of disk conformal parameterizations by over 
60% when compared with the fastest state-of-the-art approach [3]. 

2. With the significant improvement of the computational time, our proposed method possesses 
comparable accuracy as of the other state-of-the-art approaches. 

3. The bijectivity of the parameterization is supported by quasi-conformal theories. No foldings or 
overlaps exist in the parameterization results. 

4. Our proposed method is highly robust to irregular triangulations. It can handle meshes with very 
sharp and irregular triangular faces. 


4 Proposed method 

In this section, we present our proposed method for disk conformal parameterizations of simply- 
connected open surfaces in details. 

A map / : M —» N between two Riemann surfaces is called a conformal map if it preserves 
the inner product between vectors in parameter space and their images in the tangent plane of 
the surface, up to a scaling factor. More specifically, there exists a positive scalar function A such 
that f*ds 2 N = A ds 2 M . In other words, conformal maps are angle preserving. The following theorem 
guarantees the existence of several special types of conformal maps. 

Theorem 4.1 (Uniformization Theorem) Every simply connected Riemann surface is confor¬ 
mally equivalent to exactly one of the following three domains: 

(i) the Riemann sphere, 

(ii) the complex plane, 

(Hi) the open unit disk. 

Proof See [31]. □ 

With this theoretical guarantee, our goal is to efficiently and accurately compute a conformal 
map / : M —> D from a topologically disk-like surface M to the open unit disk ID). 

Before explaining our proposed algorithm in details, we point out the major differences between 
our proposed method and the two-step iterative approach [3]. Table [5] highlights the main features of 
our proposed method and the two-step iterative approach [3j for disk conformal parameterizations. 
The two-step iterative approach makes use of the disk harmonic map as an initial map, with the arc- 
length parameterized circular boundary constraint. This introduces large conformality distortions in 
the initial map. To correct the conformality distortion, two further steps are required in [3]. First, 
the Cayley transform is applied to map the initial disk onto the upper half plane for correcting the 
distortion at the inner region of the disk. Then, iterative reflections along the unit circle are applied 
for correcting the distortion near the boundary of the disk until convergence. 

In contrast, our proposed fast method primarily consists of only two stages. In the first stage, we 
find an initial planar parameterization via double covering followed by a spherical conformal map. 
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Features 

Two-step iterative 
approach [3j 

Our proposed method 

Type of input 
surfaces 

Simply-connected open 
surfaces 

Simply-connected open 
surfaces 

Initial map 

Disk harmonic map 

double covering followed 
by a fast spherical 
conformal map 

Enforcement of 
boundary when 
computing the 
initial map 

Yes 

No 

Method for 
correcting the 
conformality 
distortion 

Step 1: Use the Cayley 
transform and work on 
the upper half plane 

Step 2: Iterative 
reflections along the unit 
circle until convergence 

One-step normalization 
and composition of 
quasi-conformal maps 

Output 

Unit disk 

Unit disk 

Bijectivity 

Yes 

Yes 

Iterations required? 

Yes 

No 


Table 2 Features of the two-step iterative approach jTJl and our proposed method. 


Since there is no enforcement of the boundary in computing the initial planar parameterization, the 
conformality distortion of our initial map is much lower than that by the disk harmonic map. In the 
second stage, we enforce the circular boundary, and then alleviate the conformality distortion as well 
as achieving the bijectivity using quasi-conformal theories. The absence of iterations in our proposed 
algorithm attributes to the significant enhancement in the computational time when compared with 
the two-step iterative approach [3]. 

In the following, we explain the two stages of our proposed algorithm in details. 


4.1 Finding an initial map via double covering 


Instead of directly computing the map / : M —> B from a simply-connected open surface M to the 
unit disk B, we tackle the problem by using a simple double covering technique. The double covering 
technique was also suggested in m to compute conformal gradient field of surfaces with boundaries. 
In the following, we discuss the construction in the continuous setting. Specifically, we construct a 
genus-0 closed surface M by the following method. First, we duplicate M and change its orientation. 
Denote the new copy by M '. Then we identify the boundaries of the two surfaces: 

dM i —> dM'. (1) 


By the above identification, the two surfaces M and M' are glued along the two boundaries. Note 
that here we do not identify the interiorjaf the two surfaces M and M'. As a result, a closed surface 
is formed. Denote the new surface by M. It can be easily noted that since M and M' are simply- 
connected open surfaces, the new surface M is a genus-0 closed surface. More explicitly, denote by K 
and K g the Gaussian curvature and geodesic curvature. Assume that we slightly edit the boundary 
parts of M and M' so that M is smooth. Then by the Gauss-Bonnet theorem, we have 


/ KdA + I Kgds = 2nx(M) = 2n 
JM J dM 

[ KdA + I 

JM' J dM' 


( 2 ) 


and 


Kgds = 2ttx(M') = 2tt. 


(3) 
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Hence, we have 


2t r*(M) 



(4) 


Therefore, the new surface M has Euler characteristic y(M) = 2, which implies that it is a genus-0 
closed surface. As a remark, in the discrete case, the unsmooth part caused by the double covering 
does not cause any difficulty in our algorithm since we are only considering the angle structure of 
the glued mesh. The_details of the combinatorial argument are explained in Section [5] 

After obtaining M by the abovementioned double covering technique, we look for a conformal 
map that maps M to the unit sphere. By the uniformization theorem, every genus-0 closed surface 
is conformally equivalent to the unit sphere. Hence, the existence of such a conformal map is the¬ 
oretically guaranteed. In [2], Choi et al. proposed a fast algorithm for computing a conformal map 
between genus-0 closed surfaces and the unit sphere. The algorithm consists of two steps, and in each 
step one sparse symmetric positive definite linear system is to be solved. In the following, we briefly 
describe the mentioned spherical conformal parameterization algorithm. 

The harmonic energy functional of a map if : N —¥ § 2 from a genus-0 closed surface N to the 
unit sphere is defined as 

EW= [ IIV-0II 2 dv N . (5) 

J N 

In the space of mappings, the critical points of E(ijj) are called harmonic mappings. For genus-0 
closed surfaces, conformal maps are equivalent to harmonic maps ED- Therefore, to find a spherical 
conformal map, we can consider solving the following Laplace equation 

A T ip = 0 (6) 

subject to the spherical constraint ||^>|| = 1, where A T ij) is the tangential component of Aip on the 
tangent plane of § 2 . Note that this problem is nonlinear. In [2] , the authors linearize this problem by 
solving the equation on the complex plane: 


Acf) = 0 (7) 

given three point boundary correspondences </>(cq) = where at,bi G C for i = 1,2,3. Note that 
A T <p = Acf> = 0 since the target domain is now C. Now the problem ^ becomes linear since 
A(f> is linear and the nonlinear constraint ||'0j| = 1 in the original problem (J6| is removed. After 
solving the problem 0, the inverse stereographic projection is applied for obtaining a spherical 
parameterization. Note that in the discrete case, the conformality of the inner region on the complex 
plane is negligible but that of the outer region on the complex plane is quite large. Correspondingly, 
the conformality distortion near the North pole of the sphere is quite large. Therefore, to correct 
the conformality distortion near the North pole, the authors in [2] propose to apply the South-pole 
stereographic projection to project the sphere onto the complex plane. Unlike the result obtained by 
solving Equation 0. the part with high conformality distortion is now at the inner region on the 
plane. By fixing the outermost region and composing the map with a suitable quasi-conformal map, 
the distortion of the inner region can be corrected. Finally, by the inverse South-pole stereographic 
projection, a bijective spherical conformal parameterization with negligibly low distortions can be 
obtained. Readers are referred to j32j and (2j for more details of the harmonic map theory and the 
abovementioned algorithm respectively. 
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Fig. 1 A simply-connected open human foot model and the planar conformal parameterization obtained using 
double covering followed by spherical conformal map and stereographic projection. It can be observed that the 
resulting boundary is not perfectly circular. 


The combination of the double covering technique and the fast spherical conformal parameteriza¬ 
tion algorithm in j2] is particularly advantageous. It should be noted that because of the symmetry 
of the double covered surface, half of the entries in the coefficient matrix of the discretization of the 
Laplace Equation 0 are duplicated. Therefore, even we have doubled the size of the problem under 
the double covering technique, we can save half of the computational cost of the coefficient matrix 
by only computing half of the entries. Moreover, the spherical conformal parameterization algorithm 
in [2] involves solving only two sparse symmetric positive definite systems of equations. Therefore, 
the computation is still highly efficient. 

After finding a spherical conformal map / : M —> § 2 for the glued surface M using the param¬ 
eterization algorithm, note that by symmetry, we can separate the unit sphere into two parts, each 
of which exactly corresponds to one of M and M'. Since our goal is to find a disk conformal map 
/ : M —> ID), we put our focus on only one of the two parts. Now, we apply a Mobius transformation 
on § 2 so that the two parts become the northern and southern hemispheres of § 2 . After that, by 
applying the stereographic projection P : § 2 —> C defined by 


P{x,y,z) 




( 8 ) 


the southern hemisphere is mapped onto the open unit disk ID). Since the Mobius transformation 
and the stereographic projection are both conformal mappings, the combination of the above steps 
provides a conformal map / : M —> D. 

Theoretically, by the symmetry of the double covered surface, the boundary of the planar region 
obtained by the above stereographic projection should be a perfect unit circle. However, in the discrete 
case, due to irregular triangulations of the meshes and the conformality distortions of the map, the 
boundary is usually different from a perfect circle, as suggested in the experimental results in [3] . In 
other words, the planar region R we obtained after applying the stereographic projection may not be 
a unit disk. An illustration is given in Figure]!] To solve this issue, we need one further step to enforce 
the circular boundary, at the same time maintaining low conformality distortions and preserving the 
bijectivity of the parameterization. 
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Fig. 2 An illustration of quasi-conformal mappings. All information of a quasi-conformal map can be determined 
by the Beltrami coefficient //.. 


4.2 Enforcing the circular boundary to achieve a bijective disk conformal parameterization 


To control the conformality distortion and the bijectivity, our idea is to normalize the boundary 
and then compose the map with a quasi-conformal map. Quasi-conformal maps are a generalization 
of conformal maps, which are orientation preserving homeomorphisms between Riemann surfaces 
with bounded conformality distortions. Intuitively, a conformal mapping maps infinitesimal circles 
to infinitesimal circles, while a quasi-conformal mapping maps infinitesimal circles to infinitesimal 
ellipses with bounded eccentricity m- Mathematically, a quasi-conformal map f : C —> C satisfies 
the Beltrami equation 


df_ 

dz 


, ,df 
^ z) Tz 


(9) 


for some complex-valued functions p with ||/x||oo < 1. p is called the Beltrami coefficient of /. 

The Beltrami coefficient captures the important information of the mapping /. For instance, the 
angles and the magnitudes of both the maximal magnification and the maximal shrinkage can be 
easily determined by the Beltrami coefficient p (see Figure [2|. Specifically, the angle of the maximal 
magnification is arp(p(p))/2 with the magnifying factor 1 + |p(p)|, and the angle of the maximal 
shrinkage is the orthogonal angle (arp(p(p)) — n)/2 with the shrinking factor 1 — |p(p)|. The maximal 
dilation of / is given by: 


K(f) 


l T IIpIIoo 

1 — ||p||oo 


( 10 ) 


It is also noteworthy that / is conformal around a small neighborhood of p if and only if p(p) = 0. 
Hence, |p| is a good indicator of the angular distortions of a mapping. 

In fact, the norm of the Beltrami coefficient p is not only related to the conformality distortion 
but also the bijectivity of the associated quasi-conformal mapping, as explained by the following 
theorem: 


Theorem 4.2 If f : M —> D is a C 1 map satisfying ||/x/||oo < 1, then f is bijective. 


Proof See [TO] , 


□ 


This theorem can be explained with the aid of the Jacobian of /. The Jacobian Jf of / is given 

by 

(i - Im/I 2 )- (ii) 

Suppose \\fif\\oc < 1, then we have 1^1 0 and (1 — |p/| 2 ) > 0. Therefore, J(/) is positive 

everywhere. Since D is simply-connected and / is proper, we can conclude that / is a diffeomorphism. 
In fact, / is a universal covering map of degree 1. Therefore, / must be bijective. One important 
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consequence of Theorem |4.2| is that we can easily achieve the bijectivity of a quasi-conformal map by 
enforcing its associated Beltrami coefficient to be with supremum norm less than 1. 

Moreover, it is possible for us to reconstruct a mapping by a given Beltrami coefficient, as ex¬ 
plained by the following theorem: 

Theorem 4.3 Let Mi and M 2 be two simply-connected open surfaces. Given 2-point correspon¬ 
dences, every Beltrami coefficient p with ||/x|| < 1 is associated with a unique quasi-conformal home- 
omorphism f : Mi —* M 2 - 

Proof See m- □ 

For the aspect of numerical computations, Lui et al. [25] proposed the linear Beltrami solver 
(LBS), a fast algorithm for reconstructing a quasi-conformal map on a rectangular domain from a 
given Beltrami coefficient. The key idea of LBS is as follows. 

By expanding the Beltrami Equation |9]), we have 

= ( u * ~ Vy) + i(Vx + U y ) / 12 x 

(Ux + Vy ) 'i'ivx 

Suppose p(f) = p + ir. Then, v x and v y can be expressed as linear combinations of u x and u y : 


Vy — OL\Ux “1“ CV, 2 'U y , 

V x = Oi 2 Ux T 03'Uy, 


where 


Oil = 


012 = — 


{p-lf + p 2 
1 — p 2 — r] 2 ’ 

1 — p 2 — r) 2 ’ 

1 + 2 p + p 2 + ri 2 


«3 = 1 2 2 

1 — p z — T] z 

Similarly, we can express u x and u y as linear combinations of v x and v y : 

riy — Ol\V x H - OL c 2,V y . ) 

U x — OL2'V x + OL‘^V y . 

Hence, to solve for a quasi-conformal map, it remains to solve 

V ■ I A ( Ux U = 0 and V • (A ( Vx 1 1 = 0 


(13) 


(14) 


(15) 


(16) 


where A = 


ai OL2 
OL2 ol 3 


In the discrete case, the above elliptic PDEs (16) can be discretized into sparse linear systems. For 
details, please refer to [25] . In the following discussion, we denote the quasi-conformal map associated 
with the Beltrami coefficient p obtained by LBS by LBS(p). 

Another important property of quasi-conformal mappings is about their composition mappings. 
In fact, the Beltrami coefficient of a composition mapping can be explicitly expressed in terms of the 
Beltrami coefficients of the original mappings. 


Theorem 4.4 Let f : 17 C C —> /(I2) and g : f(f2) —> C be two quasi-conformal mappings. Then the 
Beltrami coefficient of g o f is given by 


= Pf +_ifz/fz){Pg o /) 

9 ° f 1 + ( fz/fz)pj{p g O /) 

In particular, if Pf- 1 = p g , then since Pf-i 0 / = ~{fz/fz)pf, we have 

. Pf +_(fz/fz){{-fz/lz)pf) ^ 0 

9 ° f 1 + Uz/fz)Pf(,{-fz/fz)pf) 


(17) 


(18) 


Hence g o / is conformal. 










10 


Gary Pui-Tung Choi, Lok Ming Lui 


Proof See [TO} . 


□ 


In other words, by composing two quasi-conformal maps whose Beltrami coefficients satisfy the 
above condition, one can immediately obtain a conformal map. This observation motivates the fol¬ 
lowing step. 

To enforce the circular boundary of the parameterization, we first normalize the boundary of the 
region R to the unit circle: 

'll 

(19) 


v 

V !->■ |-T 

M 


for all v G dR. Denote the normalized region by R. Since the vertices near the boundary of the region 
R may be very dense, a direct normalization of the boundary may cause overlaps of the triangulations 
as well as geometric distortions on the unit disk. To eliminate the overlaps and the distortions of R, 
we apply the linear Beltrami solver to construct another quasi-conformal map with the normalized 
boundary constraints. Then by the composition property, the composition map becomes a conformal 
map. More specifically, denote the Beltrami coefficient of the mapping g : R —y M from the normalized 
planar region to the original surface by g. We reconstruct a quasi-conformal map with Beltrami 
coefficient g on the unit disk by extending the linear Beltrami solver, so that it is applicable not only 
on rectangular domains but also circular domains. We compute a map h : R —> B by applying the 
linear Beltrami solver: 

h = LBS(/i) (20) 

with the circular boundary constraint h{v) = v for all v G dR. 

hog -1 is a conformal map from 


4.4 


Note that by the composition property stated in Theorem 
the original surface M to the unit disk B. Finally, the bijectivity of the composition map is supported 


by Theorem |4.2[ since the Beltrami coefficient of the composition map is with supremum norm less 
than 1. This completes the task of finding a bijective disk conformal parameterization. The numerical 
implementation of our proposed method is explained in Section [5} 


5 Numerical Implementation 

In this section, we describe the numerical implementation of our proposed algorithm in details. In 
the discrete case, 3D surfaces are commonly represented by triangular meshes. Discrete analogs of 
the theories on the smooth surfaces are developed on the triangulations. 

We first briefly describe the discrete version of the mentioned double covering technique for 
obtaining a genus-0 closed mesh. This discretization was also applied in nsi to compute conformal 
gradient fields of surfaces with boundaries. A triangulation K = (V, E, F ) of a smooth simply- 
connected open surface consists of the vertex set V, the edge set E and the triangular face set F. 
Each face can be represented as an ordered triple [«, v, w\ where u , v, w are three vertices. Suppose the 
boundary vertices of K are {u, ! ,}( =1 . Denote a duplicated triangulation of K by K' = (V', E', F'), 
with boundary vertices {«4}[ =1 . Let v,e,f denote the number of vertices, edges and faces of M 
respectively. We duplicate M and denote the new copy by M', with v' , e', f the number of vertices, 
edges and faces. By Euler’s formula, we have 

v — e + f = v — e + f = 1. (21) 

Then, we construct a new surface M by reversing the face orientation of M 1 and gluing the two 
boundaries dM and dM'. To reverse the orientation of K', we rearrange the order of the vertices of 
each face in F' from [m,w,w] to [u,w,v\. Then, to glue the two surfaces, we identify the boundary 
vertices of the two meshes: 

Wi <—* w'i for all z = 1,2,--- , r. (22) 

Now, denote the number of vertices, edges and faces of M by v, e and / respectively. It follows that 

v-e + f 

= (■v + v' - r) - (e + e' - r) + (/ + /') 

= v + v' — e — e' + f + f' 

= 2 , 


( 23 ) 
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Hence, the new surface M is a genus-0 closed surface. The complete double covering procedure is 
described in Algorithm [I] 


Algorithm 1: The double covering technique in our initial step. 

Input: A triangulation K = (V,E,F) of a simply-connected open surface. 
Output: A genus-0 closed mesh K = ( V,E,F ). 

1 Duplicate K and denote the copy by K' = (V ', E', F')\ 

2 Change the order of the vertices of each face in F' from [u, v, w\ to [m, w, v\\ 

3 Replace the boundary vertices w[ by Wi in E' and F' for i m 1, 2, ■ • • , r; 

4 set v = vuc'\KK=i; 

5 Set E = E U E' ; 
e Set F = E U E'- 


Then, we introduce the discretization about harmonic mappings used in the fast spherical con¬ 
formal parameterization algorithm in [2]. Let K be the triangulation of a genus-0 closed surface N. 
Let us denote the edge spanned by two vertices it, v on K by [it, v\. The discrete harmonic energy of 
ip : K —> § 2 is given by 

E W) = \ k uv\\*l>(u) - ip(v)\\ 2 , (24) 

\u,v] £K 


where k U v = cot a + cot/3 with a, /3 being the angles opposite to the edge [m, f]. This is known as 
the cotangent formula [28]. With the discrete harmonic energy, we can appropriately discretize the 
Laplace-Beltrami operator as 


Alp{Vi)= ^2 k ViVj(lp(Vj) - 1p(Vi)), 


(25) 


where N(vi) denotes the set of the 1-ring neighboring vertices of Vi. Therefore, the Laplace Equation 
0 becomes a linear system in the form 

Az = b, (26) 

where A is a square matrix satisfying 


Aij — 


k Vi vj if 7 ~ 

— J2v t eN(v i ) kv i v t if * = 3 


for all i,j. 


(27) 


For the choice of the points oi, ct 2 , « 3 , bi, 62,63 in the boundary constraint of Equation ([7]), we choose 
a triangular face [ 01 , 02 , 03 ] on the triangulation K and a triangle [ 61 , 1 ) 2 , 63 ] on the complex plane 
that shares the same angle structure as [ 01 , 02 , 03 ]. Note that the above linear system is sparse and 
symmetric positive definite. Therefore, it can be efficiently solved. 

It is noteworthy that due to the symmetry of the double covered surface, for every edge [it, v] in 
the triangulation K, there exists a unique edge [u ,v'\ in the duplicated triangulation K' such that 


k uv = cot a + cot /3 = cot a + cot ft = k u ' v >, 


(28) 


where a,/3 are the angles opposite to the edge [u, d] and a! , ft are the angles opposite to the edge 
[it , ,D / ]. Therefore, only half of the vertices and faces are needed for computing the whole coefficient 
matrix A to solve the Laplace Equation 0. More explicitly, Equation (|26[) can be expressed as the 
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following form: 



zi 


Z\ V \-r 

Cl 


Cr 

z'i 

\ Z '\V\-r) 


= b. 


(29) 


Here, {zj}^ r and T are respectively the coordinates of the non-boundary vertices of K and 

K', {Ci}i=i are the coordinates of the glued vertices, and B is a (j V| — r) x (j Vj — r) sparse symmetric 
positive dehnite matrix. It follows that we can save half of the computational cost in finding all the 
cotangent weights k uv in A. Hence, the computation of the spherical conformal map is efficient even 
the number of vertices and faces is doubled under the double covering step. 

Another important mathematical tool in our proposed algorithm is the quasi-conformal mapping. 
Quasi-conformal mappings are closely related to the Beltrami coefficients. It is important to establish 
algorithms for computing the Beltrami coefficient associated with a given quasi-conformal map, as 
well as for computing the quasi-conformal map associated with a given Beltrami coefficient. 

We first focus on the computation of the Beltrami coefficients. In the discrete case, suppose 
Ki , K '2 C R 3 are two triangular meshes with the same number of vertices, faces and edges, and 
/ : K\ —> K -2 is an orientation preserving piecewise linear homeomorphism. It is common to discretize 
the Beltrami coefficient on the triangular faces. To compute the Beltrami coefficient /i/ associated 
with /, we compute the partial derivatives on every face T\ on K\. 

Suppose Ti on K\ corresponds to a triangular face T -2 on K 2 under the mapping /. The approx¬ 
imation of /i/ on Ti can be computed using the coordinates of the six vertices of Ti and T 2 . Since 
the triangulations are piecewise linear, we can place T\ and T 2 on R 2 using suitable rotations and 
translations to simplify the computations. Hence, without loss of generality, we can assume that T\ 
and T 2 are on R". Specifically, suppose T± = [ai -H 61 , 02 +* 62 , 03 +i 63] and X 2 = [wi, u> 2 , 103], where 
ai, ( 12 , 03 , 61 , 62,63 G R, and wi,u>2,W3 G C. Recall that || = \ and || = | . 

Hence, to discretize the Beltrami coefficient, we only need to compute ^ and ^ on every triangular 
face T\. It is natural to use the differences between the vertex coordinates for the approximation. We 
define 

1 / 63 - 62 \ 1 / a 3 - a 2 

D x = — -—— 61-63 and D y = -—— 01-03 

2Area(T 1 ){ b2 _ bi J 2Area(T 1 ) ^ _ ai 



Then, we can compute the Beltrami coefficient 


M/(z) 


df /df 

dz / dz 


(31) 


on T l by 


M/CTi) 


\ ( D x + i Dy) 



I (A 


i Dy) 



(32) 


This approximation is easy to compute. Hence, it is convenient to obtain the Beltrami coefficient 
associated with a given quasi-conformal map in the discrete case. A relatively complicated task is to 
compute the quasi-conformal map / associated with a given Beltrami coefficient p./. To achieve this, 
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we apply the LBS [25] to reconstruct a quasi-conformal map from a given Beltrami coefficient, with 
the boundary vertices of the disk fixed. We now briefly explain the discretization of the LBS. 

Recall that the quasi-conformal map associated with a given Beltrami coefficient can be obtained 


by solving Equation (16). The key idea of LBS is to discretize Equation (161 into sparse SPD linear 


systems of equations so that the solution can be efficiently computed. 

For each vertex Vi , let N t be the collection of the neighboring faces attached to Vi. Let T = 
[vi,Vj,Vk] be a face ,vh = f(vt ), Wj = f(vj) and Wk = f(vk)- Suppose vi = gi+i hi and wi = si+iti, 
for l = i,j , k. We further denote the Beltrami coefficient of the face T obtained from Equation (32) 


by = pr + i r]T- Equation (14) can be discretized as follows: 


ai(T) = 


a 2 (T) = - 


(pr — 1 )~ + VT 
1 - Pt - Vt 
2t]t 


a 3 (T) = 


1 ~ Pt ~ Vt ’ 

1 + 2 pT + Pt + t] t 

1 — PJ, — IJrp 


(33) 


Then, Lui et al. [25] introduced the discrete divergence to compute the divergence operator. Define 

AT, by 


'Ai\ x (hj-hk 

Xi = Area(T 


B 


and 


1 


Arm{T) 



(34) 


Then the discrete divergence of a discrete vector field V = (Vi, V 2 ) on the triangular faces can be 
defined by 

Div(V)(vi)= J2 A *Vi(T) + bTv 2 {T). (35) 

TeNi 

With the discrete divergence, Equation (|16|) can be discretized into the following linear system: 


(36) 


Af [ai(T)or + a2(T)&T] + Bf [a.2(T)aT + a 3 (T)bx\ — 0 

TeNi 

y Aj[ai(T)c T + a2(T)d,T] + bJ[oi 2 (T)ct + a 3 (T)d,T\ = 0 

TeNi 

for all vertices u,;. Here ot, &t, Ct and dx are certain linear combinations of the ^-coordinates and 
y-coordinates of the desired quasi-conformal map /. Hence, we can obtain the x-coordinate and y- 


coordinate function of / by solving the linear system in Equation (36). For more details, please refer 
to [25]. 

With the aforementioned discrete analogs and algorithms, we are ready to describe the numerical 
implementation of our proposed algorithm for disk conformal parameterizations of simply-connected 
open surfaces. 

Given a triangular mesh M of disk topology, we first apply the double covering technique as 
described in Algorithm [T] Then, we apjrly the fast spherical conformal parameterization algorithm 
[2] on the glued genus-0 closed mesh M. To obtain a planar parameterization from the obtained 
parameter sphere S, we apply a Mobius transformation to locate the part corresponding to one of 
the two copies on the southern hemisphere. After that, by the stereographic projection, the Southern 
hemisphere is mapped to a planar region R in the complex plane. Note that this planar region 
is usually a bit different from a perfect disk due to the discretization and approximation errors. 


Therefore, we normalize the boundary of R to enforce a circular boundary as in Equation (19). 


However, this step may causes overlaps as well as conformality distortions on the normalized region 
R. 

To o vercome these issues, we apply the idea of the composition of quasi-conformal maps in Theo¬ 
rem 4.4 We first compute the Beltrami coefficient of the quasi-conformal map g : R —> M by solving 
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Equation (32) on all triangular faces of R. Denote the Beltrami coefficient by /x. We then reconstruct 
another quasi-conformal map h : R — > B with the given Beltrami coefficient /x, using LBS [25]. Note 
that boundary constraints are needed in solving Equation (36 1 . In our case, we give the Dirichlet 
condition on the whole boundary of the normalized disk R. More explicitly, the boundary condition 
used in solving Equation (36) for the map h : R — > B is 


h(v) = v 


(37) 


for all boundary vertices v G dR. After obtaining h, we can conclude that the composition map 
/ = hog -1 is the desired disk conformal parameterization. The complete implementation of our 
algorithm is described in Algorithm [2] 


Algorithm 2: Our proposed linear formulation for disk conformal parameterization 
Input: A simply-connected open mesh M. 

Output: A bijective disk conformal parameterization / : M —> B. 

1 Double cover M and obtain a genus-0 closed mesh M ; 

2 Apply the fast spherical conformal parameterization [2] on M and obtain the parameter 
sphere 5; 

3 Apply a Mobius transformation on S so that the original surface M corresponds to the 
southern hemisphere of S; 

4 Using the stereographic projection for the southern hemisphere of S, obtain a planar region R; 

5 Normalize the boundary of R: 

Vh ^w\ (38) 

for all v G dR and denote the normalized region by R; 

6 Compute the Beltrami coefficient of the map g : R —» M and denote it by p; 

7 Apply the linear Beltrami solver m to obtain a map h : R —» B 


h = LBS(/x) (39) 

with the boundary constraint h(v) = v for all v G dR. The composition map f = ho g _1 is the 
desired disk conformal parameterization; 


6 Experimental results 

In this section, we demonstrate the effectiveness of our proposed algorithm using a number of 3D 
simply-connected open meshes. The meshes are freely adapted from mesh repositories such as the 
AIM@SHAPE Shape Repository |38], the Stanford 3D Scanning Repository [39] and the Benchmark 
for 3D Mesh Segmentation !|0]- Our proposed algorithm and the two-step iterative approach [3] are 
implemented in MATLAB. The sparse linear systems in our proposed algorithm are solved using the 
backslash operator (\) in MATLAB. For the other existing algorithms in our comparison, we adopt 
available online software for the computation. All experiments are performed on a PC with a 3.40 
GHz quad core CPU and 16 GB RAM. 

We apply our proposed algorithm on various kinds of simply-connected open surfaces with dif¬ 
ferent geometry. Figure [3] shows the disk conformal parameterization of the foot model in Figure [I] 
using our proposed method. Figure [4] and [5] respectively show a simply-connected open statue model 
and a Stanford bunny model, together with the disk conformal parameterizations of them obtained 
by our proposed method. For a better visualization of the parameterization results, the triangular 
faces of the meshes are colored with the mean curvatures of the models. The disk parameterizations 









Linear Disk Conformal Parameterization 


15 



Fig. 3 


The disk conformal parameterization of the foot model in 


Figure 


q 


using our proposed algorithm. 




Fig. 4 A simply-connected open statue model and the disk conformal parameterization obtained by our proposed 
algorithm. The color represents the mean curvature of the model. 


well resemble the original models locally. The histograms of the angular distortion of our proposed 
algorithm are shown in Figure [6] It can be observed from the peaks of the histograms that the angle 
differences highly concentrate at 0. This indicates that our proposed algorithm produces disk confor¬ 
mal parameterizations with minimal distortions for different kinds of simply-connected open surfaces. 
As for the efficiency, our proposed algorithm takes only around 1 second for moderate meshes. The 
computation can also be complete within half a minute for dense meshes. 

To quantitatively assess the performance of our proposed algorithm for disk conformal parame¬ 
terizations, we compare our proposed algorithm with the existing methods. In particular, the com¬ 
parisons performed in [3] suggest that the holomorphic 1-form method m and the two-step iterative 
approach [3] achieve the best conformality among all state-of-the-art approaches as well as the bijec- 
tivity. Therefore, the holomorphic 1-form method and the two-step iterative approach are considered 
in our comparisons. The implementation of the holomorphic 1-form method is included in the Rie- 
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Fig. 5 A simply-connected open Stanford bunny model and the disk conformal parameterization obtained by our 
proposed algorithm. The color represents the mean curvature of the model. 


mannMapper Toolkit m written in C++ and the two-step iterative approach is implemented in 
MATLAB. The error thresholds in both the holomorphic 1-form method and the two-step iterative 
approach are set to be e = 10 -5 . 

The statistics of the performances of the algorithms are shown in Table[3] It is noteworthy that the 
angular distortion of our proposed method is comparable and sometimes better than the two state- 
of-the-art algorithms mm- This implies that our algorithm successfully produces low-distortion 
parameterizations. Moreover, our proposed method demonstrates a significant improvement in the 
computational time of the disk conformal parameterizations. Specifically, our proposed method is 
around 20 times faster than the holomorphic 1 -form method on average. Also, our proposed method 
requires 60% less computational time than the two-step iteration on average. The results illustrate 
the efficiency of our proposed algorithm. As a remark, in all experiments, our proposed algorithm 
generates folding-free parameterization results. 

Our proposed method is very robust to irregular triangulations. Figure [7] shows a simply-connected 
mask model which has a very irregular triangulation. Sharp and irregular triangles can be easily 
observed in a zoom-in of the model. Note that the holomorphic 1-form method m fails to compute 
the disk conformal parameterization of this model while our proposed method succeeds (please refer 
to Table [3J. This demonstrates the robustness of our proposed method. 

In addition, to illustrate the accuracy of our proposed algorithm, we map a rainbow checkerboard 
texture onto a simply-connected bunny mesh using our proposed parameterization algorithm. In 
Figure[ 8 j it can be easily observed that the orthogonal checkerboard structure is well preserved under 
the our proposed parameterization scheme. The preservation of the orthogonality indicates that our 
proposed algorithm is with negligibly low conformality distortion. 

Because of the efficiency, accuracy, bijectivity and robustness, our proposed algorithm is highly 
practical in applications, such as texture mapping in computer graphics and fashion design. For 
instance, clothing designs can be easily visualized using our proposed algorithm. Figure [9] shows a 
T-shirt model with a 2D flower pattern design mapped onto it. It is noteworthy that the patterns 
are well preserved on the 3D T-shirt model. Unlike the traditional 2D images of clothes, the 3D 
models with our texture mapping technique provide a more comprehensive and realistic preview of 
the clothes for the designers. In addition, our proposed algorithm can be used in developing virtual 
dressing rooms. In virtual dressing rooms, it is desirable to have an efficient and accurate way for 
customers to virtually try on clothes. Our proposed algorithm can help creating such a virtual dressing 
experience for online shoppers. 

As a final remark, we discuss the possibility of further boosting up the efficiency of our proposed 
algorithm. Recall that in our proposed algorithm, we need to compute the spherical conformal pa- 
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Angle difference (Degrees) 


Fig. 6 The histograms of the angular distortion of our parameterization algorithm. Here the angular distortion 
refers to the difference (in degrees) between the new angle after the parameterization and its corresponding original 
angle on the mesh. Top left: The result of statue. Top right: The result of foot. Bottom: The result of Stanford 
bunny. 


rameterization of the double covered surface using the algorithm in [ 2 ] in order to find an initial 
planar parameterization. The algorithm in [2] consists of two major steps, namely a North-pole step 
and a South-pole step. In the North-pole step, the Laplace equation 0 is solved in a triangular 
domain [&i, 62 , 63 ] on the complex plane and the inverse stereographic projection is applied. Then the 
South-pole step aims to correct the conformality distortion near the North pole of the sphere caused 
by the discretization and approximation errors. In fact, since we are only interested in half of the 
glued surface, the South-pole step may be skipped as we can take the Southern hemisphere obtained 
by the first step as our result. It may already be with acceptable conformality. 

The conformality distortion in the North-pole step in [2] is primarily caused by the choice of the 
boundary triangle [ai, 02 , 0 . 3 ]. If the chosen boundary triangle and its neighboring triangular faces are 
regular enough, then the overall distortion obtained in the North-pole step is already very low and 
only a small region around the North-pole is with relatively large distortion. In this case, the South- 
pole step is not needed to improve the distortion of the Southern region which we are interested. On 
the other hand, with the presence of the South-pole step, the overall distortion in the final spherical 
parameterization result will be negligible regardless of the choice of the boundary triangle. In short, 
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Surfaces 

No. of 
faces 

Holomorphic 
1 -form |I 2 ] 

Two-step iteration 

m 

Our proposed 
method 

Time (s) / Mea,n(\distortion\) (degrees) / SD(| distortion\) (degrees) 

Horse 

9K 

fail 

0.76 / 4.58 / 5.45 

0.18 / 4.60 / 5.49 

T-shirt 

14K 

18.64 / 1.38 / 3.26 

2.09 / 1.34 / 3.25 

0.34 / 1.35 / 3.26 

Foot 

20K 

11.73 / 1.40 / 1.24 

1.75 / 1.42 / 1.22 

0.47 / 1.42 / 1.22 

Chinese 

lion 

30K 

29.87 / 1.44 / 2.05 

2.70 / 1.42 / 2.04 

0.92 / 1.42 / 2.05 

Sophie 

40K 

28.29 / 0.36 / 0.61 

5.87 / 0.34 / 0.60 

1.31 / 0.35 / 0.60 

Bimba 

50K 

28.04 / 1.29 / 1.78 

2.22 / 1.22 / 1.74 

1.32 / 1.22 / 1.75 

Human 

face 

50K 

28.45 / 0.55 / 1.84 

4.61 / 0.53 / 1.82 

1.40 / 0.53 / 1.83 

Niccolo 

da 

Uzzano 

50K 

29.49 / 0.78 / 1.73 

7.95 / 0.75 / 1.75 

1.34 / 0.76 / 1.74 

Mask 

60K 

fail 

7.08 / 0.25 / 0.33 

1.93 / 0.25 / 0.33 

Bunny 

70K 

40.14 / 1.08 / 1.80 

4.18 / 1.08 / 1.79 

1.99 / 1.08 / 1.79 

Brain 

100K 

58.49 / 1.46 / 1.59 

6.81 / 1.46 / 1.59 

2.73 / 1.46 / 1.59 

Lion vase 

100K 

93.64 / 1.44 / 1.91 

5.34 / 1.27 / 1.75 

2.64 / 1.27 / 1.76 

Max 

Planck 

100K 

75.92 / 0.61 / 0.80 

6.58 / 0.61 / 0.80 

2.88 / 0.61 / 0.80 

Hand 

110K 

63.90 / 1.40 / 1.99 

7.51 / 1.21 / 1.31 

3.30 / 1.21 / 1.31 

Igea 

270K 

173.60 / 0.40 / 0.72 

54.63 / 0.40 / 0.71 

9.47 / 0.40 / 0.71 

Julius 

Caesar 

430K 

295.63 / 0.21 / 0.67 

65.54 / 0.20 / 0.67 

19.51 / 0.20 / 0.67 


Table 3 Performance of our proposed method and two other state-of-the-art algorithms. Here the distortion refers 
to the angular distortion of the parameterization, that is, the difference (in degrees) between the new angle and its 
corresponding original angle on the mesh. 


with a well chosen boundary triangle [ 01 , 02 , 03 ] in Equation 0 , half of the computational cost in 
computing the spherical conformal mapping can be further reduced. 

Table [4] shows the performance of the current version of our proposed method and the possible 
improved version of it without the South-pole step in [ 2 ], under a suitable choice of the boundary 
triangle [ai , a 2 ,03] in Equation 0. It can be noted that the differences in the conformality distortion 
of the two versions are mostly negligible, while the version without the South-pole step possess a 
further 40% improvement in the computational time on average. However, it should be reminded 
that the possible improved version without the South-pole step requires a suitably chosen boundary 
triangle [ 01 , 02 , 03 ] for [2], while the current version of our method is fully automatic. Hence, the 
current version of our proposed method is probably more suitable for practical applications until an 
automatic algorithm for searching for the most suitable boundary triangle [01,02,03] is developed. 


7 Conclusion and future works 

In this paper, we have proposed a linear formulation for the disk conformal parameterizations of 
simply-connected open surfaces. We begin the algorithm by obtaining an initial planar parameteri¬ 
zation via double covering and spherical conformal mapping. Note that even the size of the surface is 
doubled by double covering, the combination of the double covering technique and the spherical con¬ 
formal mapping results in an efficient computation because of the symmetry. After that, we normalize 
the boundary and compose the map with a quasi-conformal map so as to correct the conformality 
distortion and achieve the bijectivity. Our proposed formulation is entirely linear, and hence the 
computation is significantly accelerated by over 60% when compared with the fastest state-of-the-art 
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Fig. 7 A mask model with a highly irregular triangulation. Sharp triangles can be easily observed in a zoom-in of 
the model. Our proposed method can handle this kind of irregular triangulations. Left: The mask model. Right: A 
zoom-in of the model. 



Fig. 8 The Stanford bunny model with a rainbow checkerboard texture mapped onto it using our proposed param¬ 
eterization algorithm. 


approaches. At the same time, our parameterization results are of comparable quality to those pro¬ 
duced by the other state-of-the-art approaches in terms of the conformality distortions, the bijectivity 
and the robustness. Therefore, our proposed algorithm is highly practical in real applications, espe¬ 
cially for the problems for which the computational complexity is the main concern. In the future, 
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Fig. 9 A T-shirt model with a flower pattern design mapped onto it using our proposed parameterization algorithm. 


Surfaces 

No. of 
faces 

Our proposed method 

Our proposed method 
without the South-pole step 
in [2] 

Time (s) / Me&n(\distortion\) (degrees) / SD(| distortion |) (degrees) 

Horse 

9K 

0.18 / 4.60 / 5.49 

0.11 / 4.63 / 5.50 

T-shirt 

14K 

0.34 / 1.35 / 3.26 

0.22 / 1.37 / 3.27 

Foot 

20K 

0.47 / 1.42 / 1.22 

0.31 / 1.43 / 1.23 

Chinese lion 

30K 

0.92 / 1.42 / 2.05 

0.56 / 1.43 / 2.05 

Sophie 

40K 

1.31 / 0.35 / 0.60 

0.81 / 0.37 / 0.62 

Bimba 

50K 

1.32 / 1.22 / 1.75 

0.75 / 1.23 / 1.76 

Human face 

50K 

1.40 / 0.53 / 1.83 

0.83 / 0.55 / 1.84 

Niccolo da 
Uzzano 

50K 

1.34 / 0.76 / 1.74 

0.78 / 0.77 / 1.74 

Mask 

60K 

1.93 / 0.25 / 0.33 

1.22 / 0.28 / 0.37 

Bunny 

70K 

1.99 / 1.08 / 1.79 

1.14 / 1.09 / 1.79 

Brain 

100K 

2.73 / 1.46 / 1.59 

1.55 / 1.49 / 1.60 

Lion vase 

100K 

2.64 / 1.27 / 1.76 

1.45 / 1.28 / 1.76 

Max Planck 

100K 

2.88 / 0.61 / 0.80 

1.54 / 0.63 / 0.80 

Hand 

110K 

3.30 / 1.21 / 1.31 

1.74 / 1.22 / 1.31 

Igea 

270K 

9.47 / 0.40 / 0.71 

5.22 / 0.54 / 0.74 

Julius Caesar 

430K 

19.51 / 0.20 / 0.67 

12.41 / 0.21 / 0.68 


Table 4 Performance of the current version of our proposed method and a possible improved version of it without 
the South-pole step in [2], provided a suitable choice of the boundary triangle in Equation 0 - 


we plan to explore more applications, such as remeshing and registration of simply-connected open 
surfaces, based on the proposed parameterization scheme. 
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